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1 . INTRODUCTION 


" >>  Acoustically  induced  fatigue  failures  in  aircraft  structural  component 

have  resulted  in  unacceptable  maintenance  and  inspection  burdens.  In  some 
cases , sonic  fatigue  failures  have  resulted  in  major  structural  redesigns. 
Therefore,  accurate  prediction  methods  are  needed  to  determine  the  sonic 
fatigue  life  of  structures. 


Many  analytical  and  experimental  programs  to  develop  sonic  fatigue  design 
criteria,  however,  have  repeatedly  sho\flu.a.J?oor  comparison  between  analytical 
and  measured  maximum  root-mean-square '(RMS)  strain  at  high  sound  pressure 
lgyelSj.  [1,2].  Deviations  in  excess  of  100  percent  are  common.  Neglecting 
/large-deflections  in  the  analysis  has  been  Identified  as  a major  reason  for 
/ the  enormous  discrepancy  [2],  Recently,  analytical  efforts  [3,4]  with  a 
/ single-mode  approach,  and  a test  program  [5]  have  demonstrated  that  the 
/ prediction  of  panel  random  response  is  greatly  improved  by  including  the 
/ large-amplitude  effect  in  the  formulation.  Teat  results  [5]  also  show  that 
there  is  more  than  one  mode  responding.  Multiple  modes  were  also  observed 
by  White  [6]  in  experimental  studies  on  aluminum  and  carbon  fiber-reinforced 
plates  under  acoustic  loading.  White  also  showed  that  the  fundamental  mode 
responded  significantly  and  contributed  more  than  one-half  of  the  total 
mean-square  strain  response.  Therefore,  in  order  to  have  an  accurate 
determination  of  the  random  response  of  a structure,  multiple  modes  should  be 
used  in  the  analysis.  YThla  paper  presents  an  analytical  solution  technique 
for  the  large-amplitude  random  response  of  clamped  rectangular  plates 
considering  multiple  modes  in  the  analysis.  ...The  von  Krfrmin  large-deflection 
plate  equations  are  solved  by  a technique  wnlch  reduces  the  fourth  order 
nonlinear  partial  differential  equations  to  a set  of  second  order  nonlinear 
differential  equations  with  time  as  the  independent  variable.  A Fourier-type 
series  representation  of  the  out-of-plane  deflection  and  stress  function  is 
assumed.  The  compatibility  equation  is  solved  by  direct  substitution,  and 
the  equilibrium  equation  is  solved  through  the  use  of  Bubnov-Galerkin 
approach.  The  acoustic  excitation  is  assumed  to  be  Gaussian.  The 
Krylov-Bogoliubov-Caughey  equivalent  linearization  method  [7,9]  is  then  UBed 
so  that  the  derived  set  of  second  order  nonlinear  differential  equations  are 
linearized  to  an  equivalent  set  of  second  order  linear  differential  equations. 
Transformation  of  coordinates  from  the  generalized  displacements  to  the 
normal-mode  coordinates,  and  an  iterative  scheme  are  employed  to  obtain  RMS 
maximum  panel  deflection,  RMS  maximum  strain,  and  equivalent  linear  (or 
nonlinear)  frequencies  for  rectangular  plates  at  various  excitation  pressure 
spectral  densities.  Convergence  of  the  resultB  is  demonstrated  by  using  4, 

6,  10,  and  15  terms  in  the  transverse  deflection  function. 

2.  MATHEMATICAL  FORMULATION  AND  SOLUTION  PROCEDURE 

Assuming  that  the  effects  of  both  inplane  and  rotatory  inertia  forces 
can  be  neglected,  the  dynamic  von  Karmen  equations  of  a rectangular  isotropic 
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plate  are 

L(w,F)  ■ DV4w  + phw,tt  + gw,t  - p(t) 

h(F,yyW,xx  + F.XXW,yy  “ ^ F . ^ ) “ 0 

74F  - E(w?w  - «,xxw.yy) 

The  transverse  deflection  which  satisfies  the  clamped  boundary 
conditions 
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is  assumed  as 

V<y,t>  ' ” 1 l 

m n 

m,n  ■ 1 ,2,3,... 

(A) 

where 

fm(x)  ■ cos  - cos 

(m+l)nx 

a 
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9n(y)  - cos  lEjfry  - cos 

(n+l)«y 
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Upon  examination  of  the  foregoing  expression  for  the  transverse 
deflection,  It  Is  found  that  the  compatibility  equation  (2)  can  be 
Identically  satisfied  If  the  stress  function  F 1b  taken  In  the  following 
form: 


Direct  substitution  of  Eqs.  (A)  and  (7)  into  Eq.  (2),  performing  the 

required  differentiations,  multiplications,  and  a Fourier  analysis  of  the 

resulting  terms  yields  a quadratic  relationship  between  F and  W 
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In  which  B , , are 
pqmnkl 


Integers  (tabulated  In  Reference  10)  and  a»a/b.  A 


complete  description  of  the  solution  technique  used  to  solve  equation  (2)  Is 
given  in  Reference  [10]. 


The  average  edge  loads  P and  P in  Eq.  (7)  are  determined  from  inplane 

x y 

boundary  conditions.  The  particular  inplane  boundary  condition  of  most 
interest  in  the  study  of  sonic  fatigue  of  structural  panels  is  the  one  in 
which  the  edges  are  restrained  from  movement,  that  is 
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v“0  on  y«0  and  b 

or 
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Performing  the  differentiation  and  integration  as  indicated  in  Eq.  (10) 
yields  relationships  for  P and  P in  terms  of  the. deflection  coefficients 

w x y 

mn'  These  relations  are  too  lengthy  to  reproduce  here,  but  may  be  found 
in  References  [101  end  [11]. 

With  the  assumed  deflection  w given  by  Eq.  (4)  and  the  stress  function  F 
given  by  Eq.  (7),  Eq.  (1)  is  then  satisfied  by  applying  Bubnov-Galerkin 
method: 

//•  L ( w,F ) f rg sdx4y  ■ 0 r,s  ■ 1,2,3,...  (11) 


The  integration  of  each  of  the  terms  in  Eq.  (11)  can  be  found  in  Ref.  [10]. 

A set  of  nonlinear  time-differential  equations  is  obtained  after  performing 
the  integration  over  the  total  area  of  the  panel,  and  it  can  be  written  in 
matrix  notation  as 

[M]{W>  + [C]{ft>  + [K]l{W>  + (P(W)>  • <p(t)>  (12) 

where  [M] , [C] , and  [K]^  are  the  generalized  mass,  damping,  and  linear 

stiffness  matrices,  respectively,  and(B}  is  a vector  function  which  is  cubic 
in  the  generalized  displacements  {W  }. 

If  the  acoustic  pressure  excitation  p(t)  is  stationary  Gaussian, 
ergodic,  and  has  a zero  mean,  then  application  of  the  Krylov-Bogoliubov- 
Caughey  equivalent  linearization  technique  [3,4, 7-9]  yields  an  equivalent  set 
of  linear  equations  as 

[M]{ii>  + [C]{5l}  + ([K]l  + [K]el)  {W}  - <p(t>)  <13a) 

or 

Cm] { w>  + [c3{fl>  + DOW  - {p( t>>  (13b) 

The  elements  of  the  generalized  equivalent  linear  or  nonlinear  stiffness 
matrix  [K]^  Ran  be  derived  from  the  expression  [8] 


<KEL5rs1J 


(14) 


where^O  Is  an  expected  value  operator.  The  elements  (Kj,^  8;  . are  too 

lengthy  to  reproduce  here,  but  may  be  found  in  Reference  [11].  The 
approximate  generalised  displacements  {W },  computed  from  the  linearized  Eq. 
(13) , are  also  Gaussian  and  approach  stationary  because  the  panel  motion  is 
stable. 


2 

To  determine  the  mean-square  generalized  displacements  W 

mn 

an  iterative  process  is  Introduced.  The  undamped  linear  equation  of  Eq. 
(13a)  is  solved  first,  which  required  simply  the  determination  of  the 
eigenvalues  and  eigenvectors  of  the  undamped  linear  equation 


wj  [M3U>j  - [K]lU>j 


(15) 


where  ui.  is  the  linear  frequency  of  vibration,  and  {$},  is  the  normal  mode 
shape.  3 3 

Apply  a coordinate  transformation,  from  the  generalized  displacements  to 
the  normal  coordinates  (this  analysis  will  use  the  first  4 modes),  by 

00  - U]  (q)  4<m  (16) 

mxl  mx4  4x1 

where  each  column  of  [<j>]  la  a normal  mode  {<$>}..  The  damped  linear  aquation  of 
Eq.  (13a)  becomes  3 

CM]{q>  + CC](q}  ♦ CKJL(q)  - (P(t)> 


where 


CM]  - U3t[m][*3 

CC]  • U3TtC3U3  - 2C*Cu)3CM] 
CKJl  - [43tCK3lC43  • Ch>2JCM3 

<P>  • C43t(p> 

The  jth  row  of  Eq.  (17)  la 

qj  + 2Cj  U.J  qj  * ^ dj  ■ ^ 

The  mean-square  normal  coordinate  is  simply 
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where  S (w)  is  the  spectral  density  function  of  the  excitation  P,(t).  The 

* J 

convariance  matrix  of  the  linear  generalized  displacements  is 
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This  Initial  estimate  of  expected  value  on  generalised  displacements  can  now 
be  used  to  compute  the  generalized  linear  stiffness  matrix  [K]_.  through  Eq. 

CtL 

(14).  The  undamped  linearized  equation  of  Eq.  (13)  Is  solved  again 

QjDOMj  - (Ck]l  + [K]EL){»>j  (22) 

where  j iB  the  equivalent  linear  or  nonlinear  frequency  of  vibration,  and 
{ } Is  the  associated  equivalent  linear  normal  mode  shape.  Then  Eq.  (13) 
is  transformed  again  and  has  the  form  as 

CttUq)  + CCj{q>  + ftO{q>  « {P(t)> 

In  which 

CKO  - U]t(CK]l  + EK]el)U]  - Cq2OCmO 

the  Jth  row  of  Eq.  (23)  is 

’j  + VjV  “I’j  'HJ 

and  the  displacement  convariance  matrix  Is 

The  diagonal  terms  of  [W^W^]  are  the  mean-square  generalised  displacements 
7mn.  As  the  Iterative  process  converges  on  the  lter-th  cycle,  the  relations 
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" ^J^lter  - 1 
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^QjNter  " ^DjNter-l 

(27b) 

wmn  Iter  ,vWlter-l 

(27c) 

become  satisfied.  Convergence  is  considered  achieved  when  the  difference  of 
the  RMS  displacements  satisfies  the  requirement 


|(*MS  w«,n>1*r  - “mnW-1  |„„-3 


<10  for  all  m,  n 


(28) 


Once  the  RMS  generalized  displacements  are  determined,  the  RMS 
deflection  and  the  maximum  RMS  strain  can  be  determined  from  Eqs.  (4),  (7) 
and  the  strain-displacement  relations  as  follows. 


For  a rectangular  plate  (a>b)  with  clamped  support  along  all  four 
edges,  the  maximum  bending  strain  occurs  at  the  extreme-fiber  perpendicular 
to  its  long  edge  at  the  midpoint. 


or 
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(29a) 


EEfm(  x>  [(n-l)2cosi-n^y  - (n+1?  co^-^1ty  ]Wmn 


(29b) 


The  membrane  strain  is  given  by 


^cy^m  " TT  (F,xx  “ vF,yy^ 


or 


e..b‘  2 


h * m ® mn  mny  mn  mnklpq  pq  ppmnk1  mn  kl 


in  which 


From 


C _ - co&L 

M (p2/«  ♦ dz«)z  * * 


V 'Vb*  V» 


the  maximum  mean-square  strain  becomes 


For  Gaussian  random  processes  with  zero-mean  we  have 

«<M1jwklwnm>  “ 0 


(30a) 

(30b) 
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^w1jwklwmnwrs^  " ^wljwk1  * ^wmnwrs^ 

+ «(W1jWmn>  + 1^<w1jwrs)  *(WklMmn> 


(35) 


and  the  maximum  RMS  strain  can  be  determined  from  Eq.  (33). 


3.  RESULTS  AMD  DISCUSSION 


Using  the  present  formulation,  the  nonlinear  response  of  square  and 
rectangular  (a*  2)  plates  with  all  edges  clamped  and  subjected  to  broadband 
random  excitation  are  studied.  In  the  results  presented,  the  white  noise 
excitation  is  band-limited  with  a frequency  bandwidth  of  25  Hz  to  6000  Hz, 
the  damping  ratio  is  assumed  to  be  constant  for  all  four  normal  modeB  and 
Poisson* s ratio  is  equal  to  0.3.  Mean-square  center  deflection  and  maximum 
mean-square  strain  are  presented  in  a nondimensional  form.  The 
nondimenaional  forcing  spectral  density  parameter  is  defined  as 
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2 TT  Sp  (to) 
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Also,  since  the  loading  is  symmetric,  only  symmetric  generalized 
displacements  W are  retained  in  the  transverse  deflection  function. 

Dul 


The  convergence  of  the  solution  technique  waa  examined  in  order  to 
determine  the  degree  of  accuracy  possible  with  a highly  truncated  transverse 
deflection  function  series.  The  mean-square  center  deflection  versus  the 
non-dimensional  spectral  density  parameter  using  4,  6,  10,  and  15  terms  in 
the  deflection  function  for  a square  plate  la  shown  in  Figure  1.  The 
particular  generalized  displacements  that  make  up  the  various  orders  of  the 
deflection  function  are  shown  in  Table  1.  Figure  1 clearly  indicates  that  a 
6-term  solution  gives  accurate  results  for  the  nonlinear  maximum  defection 
while  a 4-term  solution  will  provide  accurate  linear  results.  The  maximum 
strain  occurs  at  the  extreme  fiber  of  the  panel  and  at  the  midpoint  of  the 
long  edge.  The  direction  1b  perpendicular  to  the  edge.  Figure  2 shows  the 
maximum  mean-square  strain  versus  S.  for  the  square  plate  using  4,  6,  10  and 
15  terms  in  the  deflection  function.  The  convergence  of  the  mean-square 
strain  is  much  slower  an  compared  with  that  of  the  mean-square  deflection. 

TABLE  I 

Generalized  Displacements  Used  in  Convergence  Studies 

Number  Generalized  Displacements 

of  terms  Wn  W13  W3J  W33  W.j  V?51  W35  WJ3  W1?  W?1  WJ5  W3?  W?3  W19  W91 
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Figures  3 and  4 show  the  maximum  mean-square  nondimensional  deflection 
versus  the  nondimensional  spectral  density  of  acoustic  pressure  excitation 
for  rectangular  panels  of  aspect  ratios  of  1 and  2 with  the  damping  ratio 
equal  to  0.009,  0.018  and  0.027.  Figures  5 and  6 show  the  maximum 
nondimensional  mean-square  strain  versus  the  nondimensional  spectral  density 
for  rectangular  panels  of  aspect  ratios  of  1 and  2 with  the  damping  ratio 
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equal  to  0.009,  0.018  and  0.027.  Ten  terns  were  Included  In  the  deflection 
function,  to  generate  the  results  shown  in  Figures  3 through  6, 

4.  CONCLUDING  REMARKS 

An  analytical  solution  technique  is  presented  for  determining  the 
large-amplitude  random  response  of  clamped  rectangular  panels  while  including 
multiple  modes  in  the  analysis.  Accurate  mean-square  deflections  can  be 
obtained  with  the  use  of  6 terms  in  the  deflection  function,  while  it  is 
necessary  to  consider  as  many  as  10  or  more  terms  for  the  accurate  deter- 
mination of  the  strains.  In  the  numerical  examples  presented  , a constant 
damping  ratio  for  all  four  modes  has  been  used.  However,  nonlinear  damping 
phenomena  have  been  observed  in  experiments  [5,6].  Therefore,  further 
effort,  is  needed  to  better  understand  the  effects  of  nonlinear  damping  on 
panel  response. 
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Figure  1.  Convergence  of  the  Mean- 
Square  Center  Deflection 
for  the  Square  Plate 
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Figure  2.  Convergence  of  the  Maximum 
Mean- Square  Strain  for  the 
Square  Plate  ( £ a 0.009) 
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Figure  3.  Mean-Square  Center  Deflection 
vs  Pressure  Spectral  Density 
for  a Clamped  Square  Plate 


Figure  4.  Mean-Square  Center  Deflection 
vs  Pressure  Spectral  Density 
for  a Clamped  Rectangular 
( “■  2)  Plate 


Figure  5.  Maximum  Mean-Square  Strain 
vs  Pressure  Spectral  Density 
for  a Clamped  Square  Plate 
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Figure  6.  Maximum  Mean-Square  Strain 
vs  Pressure  Spectral  Density 
for  a Clamped  Rectangular 
( «■  2)  Plate 


10.  RESPONSE  TO  SEISMIC 
EXCITATION 


